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ABSTRACT 

The phenomena customly called Dark Matter or Modified Newtonian Dynamics (MOND) have 
been argued by Bekenstein (2004) to be the consequences of a covariant scalar field, controlled by 
a free function (related to the MOND interpolating function fL{g/af)) in its Lagrangian density. 

In the context of this relativistic MOND theory (TeVeS), we examine critically the interpolating 
function in the transition zone between weak and strong gravity. Bekenstein’s toy model produces 
too gradually varying fi and fits rotation curves less well than the standard MOND interpolating 
function fl(x) = + However, the latter varies too sharply and implies an implausible 

external field effect (EFE). These constraints on opposite sides have not yet excluded TeVeS, but 
made the zone of acceptable interpolating functions narrower. An acceptable “toy” Lagrangian 
density function with simple analytical properties is singled out for future studies of TeVeS in 
galaxies. We also suggest how to extend the model to solar system dynamics and cosmology, and 
compare with strong lensing data (see also astro-ph/0509590). 

Subject headings: gravitation - dark matter - galaxy kinematics and dynamics 

where the interpolating function fi{x) is a func¬ 
tion which runs smoothly from fL{x) = cc at a; 1 
to il{x) = 1 at a: ^ 1 with a dividing gravity 
scale oo ~ 10“®cms“^ ~ cHo/6 at the transi¬ 
tion. This simple correlation was taken as the ba¬ 
sis for the MOND theory (or more precisely the 
aquadratic Lagrangian theory) by Bekenstein & 
Milgrom (1984, hereafter BM84), where one mod¬ 
ifies the Newtonian gravity of a baryonic galaxy 
to eliminate the need for dark matter. 

Recently interests on the subject of MOND 
have been further stimulated since Bekenstein 
(2004, hereafter B04) provided a Lorentz-covariant 
theory (dubbed TeVeS), which passes standard 
tests to check General Relativity, and allows for 
rigourous modeling of Hubble expansion and grav- 


1. Introduction 

On galaxy scales, dark matter generally domi¬ 
nates over baryons (stars plus gas) at large radii. 
At intermediate radii in a galaxy where dark mat¬ 
ter and baryons overlap with comparable amounts, 
the two mass profiles are not uncorrelated (Mc- 
Gaugh 2005). The correlation between the New¬ 
tonian gravity of the baryons and the over¬ 
all gravity g (baryons plus dark matter) can be 
loosely parameterized by the Milgrom’s (1983) em¬ 
pirical relation 

A(s/ao)g = gN, (1) 
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itational leasing (e.g. Zhao et al. 2005). In 
TeVeS the MONDian behaviour originates from 
a scalar field, the dynamics of which is controlled 
by a Lagrangian density involving a free func¬ 
tion that yields the expected dynamics in the 
low-acceleration limit (although BM84 theory is 
not precisely recovered). This freedom of the La¬ 
grangian density, that echoes the freedom in the 
choice of the interpolating function jl in MOND, 
means that every choice of the free function defines 
a distinct theory. As this class of theories do not 
at present derive from any basic principle and are 
purely phenomenological, the only constraints on 
the free function must come from phenomenologi¬ 
cal grounds. A refinement of the function studied 
by B04 as a toy model is surely needed. In this let¬ 
ter, we differentiate popular choices of the MOND 
pL function by fitting a benchmark rotation curve, 
and argue that many of those functions are likely 
unphysical in the TeVeS context. We then pro¬ 
pose a new free function for TeVeS in the domain 
relevant for galaxies, with a possible extension to 
solar system dynamics and cosmology. 

1.1. Warming up to TeVeS 

It is a tensor-vector-scalar Lorentz-covariant 
field theory, where the tensor is the Einstein met¬ 
ric gafi out of which is built the usual Einstein- 
Hilbert action, Ua is a dynamical normalized vec¬ 
tor field, and (j) a dynamical scalar field. The ac¬ 
tion is the sum of the Einstein-Hilbert action for 
the tensor g^fj, the matter action, the action of the 
vector field Ua, and the action of the scalar field 
(f>. Einstein-like equations are obtained for each 
of these fields by varying the action w.r.t. each of 
them. The action of the scalar field (p involves a di¬ 
mensionless parameter k (of the order of a few per¬ 
cents), a length scale parameter I (~ 
a free dimensionless function linking oc y 

with the auxiliary nondynamical scalar field g. 

More relevant to us, the physical metric in 
TeVeS near a quasi-static galaxy or the solar sys¬ 
tem is identical to that of General Relativity, with 
a potential 

$ = S$Ar -k (p (2) 

where 5 ~ 1. This means that the scalar field 
(p plays the role of the dark matter gravitational 
potential. It is related to the Newtonian potential 
4 )N (generated by the baryonic density p) through 


the equation (similar to the field equation for the 
full $ in BM84) 

=47rGp, (3) 


where ps is a function of the scalar field strength 
gs = |V^|. It is related the p function of B04 and 
the interpolating function p of MOND by 



-M, t 's Jl_ 

1 — p 47r 


( 4 ) 


In the intermediate to deep-MOND regime, the 
toy free function in the scalar action of B04 gives 
rise to the following interpolating function: 


p{x) 


VI + 4a;- 1 

VTT^+i 


( 5 ) 


2. Rotation Curves 


Milgrom’s formula (1) is a good approxima¬ 
tion of the BM84 and TeVeS theories since for 
most disk galaxies the curl field is negligible or ex¬ 
actly zero when solving Eq. (3) (Brada & Milgrom 
1995). It provides a plausible unified picture for 
the phenomenology of individual galaxies, which is 
more challenging to understand in the dark mat¬ 
ter framework. However, while the precise func¬ 
tional form of the interpolating function is not nec¬ 
essary to make many fundamental predictions (see 
Sanders & McGaugh 2002), it is nevertheless im¬ 
portant in order to fit the rotation curves of galax¬ 
ies. Indeed one of the most striking successes of 
MOND is the ability of Milgrom’s formula to fit 
the rotation curves of a wide range of galaxies with 
virtually the same value of ao = 1.2 x 10“®cms“^, 
and the same “standard” interpolating function: 


p(x) 


X 

vmv 


( 6 ) 


This “standard” interpolating function was origi¬ 
nally put in by hand, and does not derive from any 
physical principle. However, with the number of 
galaxies with good data increasing, the freedom of 
this function should be restrained by the observa¬ 
tions. Famaey & Binney (2005, hereafter FB05) 
have e.g. found that the “simple” interpolating 
function. 


p{x) 


X 

1 + x’ 


( 7 ) 


gives a better fit to the terminal velocity curve of 
the Milky Way than Eq. (6), while yielding an 
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extremely good fit to the rotation curve of the 
standard external galaxy NGC3198 (cf. Fig. 1). 
In comparison the toy model in B04 gives rise to 
Eq. (5) in spherical symmetry in the intermedi¬ 
ate to deep-MOND regime. This interpolating 
function triggers too slow a transition from the 
MONDian to the Newtonian regime in the bench¬ 
mark rotation curve of NGC3198 (cf. Fig. 1); the 
same conclusion is reached for the terminal veloc¬ 
ity curve of the Milky Way (FB05). 

In short, from the analysis of rotation curves, 
Eq. (6) and Eq. (7) are preferred over Eq. (5). 
FB05 noted that Eq. (7) is preferred over Eq. (6) 
to fit the TVG of the Milky Way. They also de¬ 
rived the best MONDian fit (in the strong and 
intermediate regimes only) of the Milky Way, 
and found the interpolating function to transition 
smoothly from Eq. (7) at a; < I to Eq. (6) at 
a; > 10. Of course, it is not obvious that all the 
other galaxies will give the same answer as the 
Milky Way and NGG3198. The observational con¬ 
straints presented in this section should thus be 
considered as an indication more than a rigorous 
constraint. 


3. The external field effect 


In TeVeS the potential consists of two parts 
(Eq. 2): the Newtonian potential, and the scalar 
field which satisfies the BM84 formulation (Eq. 3). 
One of the key features of the BM84 theory is the 
existence of the dilation effect of an external field, 
which is why MOND does not satisfy the strong 
equivalence principle. Gonsider the perturbation 
generated by a low-mass m body inside a domi¬ 
nating uniform external scalar field strength g^z, 
Eq. (3) for (j) becomes 


V.[^5V^] « (Aid^ + dl + d‘^)(l) = ATrGmS{r), 

( 8 ) 

which differs from the Newtonian Poisson equation 
by the dilation factor Ai given by (BM84; Zhao 
2005; Zhao & Tian 2005) 


Ai = l + 


d\ngf 


Ms = (9) 


As in BM84, the solution is an equal potential 
surface given by 


<l>{x,y,z) 


-9sZ- 


Gm 


f = z"^ + Ai ( 2/2 + a;2), 

( 10 ) 


where m /is the effective mass of the satellite 
and r is the effective distance from the centre of 
the satellite. Hence the perturbation scalar field 
is stretched in the z direction by a factor I/-/A 1 . 
To exclude an imaginary dilation, a theory should 
have 


Ai 


dln/Xsgs 
din gs 


( 11 ) 


This requires models where gs is an increasing 
function of gy = flags- The factor Ai also deter¬ 
mines the shape of the Roche Lobe of a rotating 
satellite, which have a middle-to-long axis ratio 


y 3 ^ (Zhao 2005; Zhao & Tian 2005). In fact the 

stretching factor enters almost all processes 

involving satellites (e.g. Brada & Milgrom 2000). 
For these reasons we consider models with a neg¬ 
ative Ai unphysical. In the original BM84 theory, 
this is not a problem because I < Ai < 2 thanks 
to a monotonically increasing fl{g) (not fisigs) as 
we are concerned with), so the external field effect 
is a mild curiosity. This, however, is not the case if 
the MOND effect is created by a scalar field. From 
Eq. (4), we have that models with the standard p, 
function (Eq. 6) yield a gs which increases with g^ 
to some point, and then starts decreasing in the in¬ 
termediate regime (see Fig. 2). At the same scalar 
field strength gs, there are two different (spherical) 
Newtonian gravity strength, i.e., the scalar func¬ 
tion fisigs) becomes multi-valued for the same gs, 
hence is ill-defined.^ This is a general feature of 
any sharply increasing MOND p function (e.g., the 
“exponential” function p = 1 — exp(—cc), and the 
best fit of FB05). This type of behaviour is unde¬ 
sirable in a physical model for the scalar field. On 
the other hand the simple function (Eq. 7) and 
B04 toy model (Eq. 5) give physical monotonic 
fisigs), hence positive Ai. 


4. New free function 

As a toy model B04 chose for the /x-function (see 
Eq. (4) above) the implicit, discontinuous formu¬ 
lation: 

yjfl) _ ^ fi^ ifi- Af 

3fc'2 - al fc'2 4(1 - /x) ’ ^ ’ 

^The ^s{s) is given by two root branchs of an essentially 
4-th order polynomial equation s = 
a very lengthy multi-valued function. 
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where notations are as in B04. Here cosmology 
> 2) and galaxies (0 < ^ < 1) are completely 
detached from each other, while the fit to galactic 
rotation curves is poor. 

Our aim here is to propose a new explicit, 
monotonic and continuous interpolating function 
A simple choice could be 



_M_ 


n = 0,1. 


(13) 

Here y > 0 and 0 < /i < 1 corresponds to quasi¬ 
static systems as in B04, and y < 0 to cosmol¬ 
ogy. This toy function is consistent with rotation 
curves, since it reduces to the simple interpolat¬ 
ing function Eq. (7) in the range for galaxies, i.e., 
p, = k'ps ~ (0 — 10)fc' <C 1, where we have the 
explicit relations 


Ps 




(14) 


which are ready to be fed into a solver for Eq. (3). 


4.1. Solar system 


Our toy function also recovers solar system dy¬ 
namics at a level similar to B04, if not better. 
For example, in n = 1 case (the default case), 
the scalar field goes to infinity when p approaches 
one. This new “toy” function is easier to use 
than B04 because it allows for an analytic for¬ 
mulation for the scalar field interpolating function 
Ps{s) = -p- = in both strong and weak gravity, 
namely 


k'ps{s) 


1-k' 

2 


-1 +^[{l + k')s + l]^-4s 


2s 


(15) 


The dilation factor Ai = 1 -|- ~ 

(1 — 4) is in the plausible range for k' = 0.03 — 1. 
As we can see from Fig. 3, the new function yields 
a monotically increasing scalar field strength 
and the correction for the solar system dynamics 
is less than in B04 toy model (when adopting k = 
0.01, the variation from Sun to Saturn is an order 
smaller), which may or may not be relevant to the 
Pioneer Anomaly. 


4.2. Cosmology 

Although the B04 discontinuous cosmological 
branch of the function y{p) could as well be 


appended to our function in the range p > 2, 
it is more desirable to be in a universe where 
the Lagrangian density has a smooth transi¬ 
tion between weakly gravitating quasi-static sys¬ 
tems {y (X (V(/))^/ag > 0) and cosmology {y oc 
—2{dt4>)‘^ < 0). A possible way, as done here 
(Eq. 13), is to copy our function in y > 0 regime 
into the y < 0 regime by a simple mirror-imaging, 
(cf. Fig. 4). This means that the outer parts 
of galaxies would connect smoothly into the cos¬ 
mological expansion as y passes from 0-1- to 0—. 
While a mirror-image with fine tuning is attrac¬ 
tive, it is not a necessary condition: any posi¬ 
tive continuous function p{y) going through, e.g., 
y(0—) = p{0+) ~ 0 is worth exploring. The 
next question is whether this kind of cosmology 
can produce realistic Hubble expansion and the 
Cosmic Microwave Background. 

5. Conclusion 

Part of the amazing successes of the non- 
relativistic version of MOND in explaining galaxy 
dynamics is due to its “standard” interpolating 
function p = xj\/l + x^. When exploring a range 
of other empirical functions in the context of 
TeVeS, we see that the fit to rotation curves be¬ 
comes poorer for more gradual p function, whilst 
an external field effect with imaginary dilation 
happens for more rapid changes of the p function. 
These two independent constraints from opposite 
sides suggest a fairly narrow range of TeVeS free 
functions. Among these we propose a simple ex¬ 
pression which works for both very weak and very 
strong gravity (Eq. 15), with a possible extension 
to cosmology. Unlike the one in B04, the new func¬ 
tion has the nice feature that it links quantities 
of TeVeS and of MOND easily, hence facilitates 
future examinations using galaxy dynamics and 
solar system data. The explicit simple monotonic 
interpolating function Ps{gs/o,o) could be easily 
fed into a numerical solver for Eq. (3) (e.g., as 
developed by Ciotti et al. 2005), and could allow 
for the modelling of realistic galaxy geometries 
(the curl-field of BM84, neglected here following 
conventional wisdom, could be put back with a 
realistic amount). Combined with a galaxy with 
sensitive kinematical data, this may confirm or fal¬ 
sify our “toy” function, hence further establish or 
squeeze the parameter space of the TeVeS theory. 
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As two final remarks, (i) we note that multiple- 
imaged gravitational systems present a challenge 
to all MOND/TeVeS interpolating functions (cf. 
Fig. 5). Zhao et al. (2005) found that ellipti¬ 
cal galaxies of comparable luminosity and redshift 
show a large scatter in their Einstein ring sizes. 
Among the previously proposed ^ functions the 
B04 toy function is the most effective in lensing, 
but none of the functions seems to fit all lenses 
in the point lens case (shaded zone); the fit is 
poorer with realistic lens mass profile (cf. Fig. 
17 of Zhao et al. 2005). (ii) We also note that the 
dark matter potential is fundamentally different 
from the scalar field; although the two are some¬ 
times degenerate in fitting rotation curves, there 
is no equivalent of EFE in dark matter, hence the 
dark matter potential enjoys more freedom. 
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Qs/ao > g/ao 


Fig. 2.— shows X = g/ao vs. il{x) (lower right) 
for popular choices: “standard” (Eq. 7, thick 
dashed), “exponential” (dashed-dotted), “Milky 
Way” (fit of FB05, dotted) and “simple” (Eq. 8, 
solid) The corresponding scalar field s = (/s/oo 
vs. its modification function fXsis) is also shown 
(upper left). Circles mark where s « 0.25. Note 
popular il{x) functions often lead to multi-valued 
fis{s) curves except for “simple”. 
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Fig. 3.— Compares our proposed (solid, marked 
ZF) model with the B04 model (dashed) in the pa¬ 
rameter space of the rescaled scalar field strength 
s = gs/ao vs. Newtonian field strength gN/ao = 
sgs{s) in cases of fc = drrfc' = 1 (thick) and 
k = 0.01 (thin). Also overplotted is the deviation 
from the oc force in the solar system (flatter 
dotted line for our model, other dotted line for 
B04). This should be compared with constraints 
from planets and the measurement of the Pioneer 
Anomaly (long and short horizontal thick lines). 







Fig. 4.— The upper panel presents our newly pro¬ 
posed function /Xs(s) (Eq.(16), solid line), com¬ 
pared to B04 toy function (dashed). We also 
present a possible extension to the cosmological 
regime y < 0 by mirror-imaging (this kind of 
extension is far from being unique). The lower 
panel displays the above extention of our new La- 
grangian density L from the weak gravity y = 0-1- 
into cosmology y < 0. B04 Lagrangian is shown 
as dashed line for comparison. Note y oc (V^)^/aQ 
is large and positive for the solar system, and 
y (X is negative for cosmology. 







Fig. 5.— Error bars show the stellar Newtonian 
gravity vs. the light-bending gravity for CAS- 
TELS sample of galaxy-quasar lenses (Zhao et al. 
2005), where we allow for a factor of two uncer¬ 
tainty with the stellar mass-to-light ratio. Pre¬ 
dictions are made for point lenses with a /i func¬ 
tion in between the B04 toy model (lower bound¬ 
ary of the hatched zone) and a fx with a sudden 
transition (upper boundary); models in the upper 
hatched zone imply unphysical scalar field. Galax¬ 
ies would require more mass to bend the light than 
the point-lens predictions here. 
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